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Abstract 

We present an exact calculation of the Luttinger liquid relation for the 

one-dimensional, two-component SC model in the interaction strength range 

— 1 < s < by appropriately varying the limits of the integral Bethe Ansatz 

equations. The result is confirmed by numerical and conformal methods. By 

a related study of the transport properties of the SC model, we can give an 

exact formula for the susceptibility. Our method is applicable to a wide range 

of models such as, e.g., the Heisenberg-Ising chain. 
72.15Nj,05.30.-d,75.30Ds 



1 



Typeset using REVTpjX 



The Luttinger liquid approach U to one-dimensional quantum models in conjunction 
with methods of conformal field theory is extremely useful for understanding the low- 
energy physics of these models. The essence of the method is the realization that close to 
the fermi level, all excitations are essentially hydrodynamic density fluctuation plus finite- 
size terms conveniently described by conformal finite size formulas ||. Furthermore, the 
finite-size spectrum may be given in terms of a single parameter due to the presence of 
the so-called Luttinger liquid relation. In most models, this parameter is the renormalized 
interaction strength. 

Although the set of models that fall into the Luttinger liquid universality class (the 
c = 1 universality class of the Gaussian model of conformal field theory) incorporates all 
"classic" short-ranged ID models such as Heisenberg-Ising and Hubbard model there is 
no rigorous proof that all possible ID models should, too. Indeed, the calculation of the 
central charge c in long-ranged models of the 1/r 2 type yields contradictory results |J . 

In this Letter, we will give an exact derivation of a Luttinger liquid type relation in a 
two-component quantum system with hyberbolic interaction defined by the Hamiltonian 

1 d 2 

ff = ~9 E o3+ E Vjkfe-Xk), (1) 

l<j<N UX j l<j<k<N 

with pair potential given as 

1 + (TjCTk 1 — CfCfc 



v jk (x) = s(s + 1) 



2sinh 2 fx) 2cosh 2 (x) 



s > -1. (2) 



The quantum number a = ±1 distinguishes the two kinds of particles. We may usefully 
think of it as either representing charge or spin. The model has been solved in Ref. || by an 
asymptotic Bethe Ansatz calculation and for large distances, its correlation functions have 
been calculated J7|. We call it the SC model for the sinh-cosh interaction. The SC model is 
closely related to the Hubbard and the Heisenberg-Ising (H-I) model || and may be viewed 
as a H-I fluid, with mobile spins . 

We restrict our investigation to — 1 < s < 0, where the system has two gap less excitations 
corresponding to a particle-hole and a two spin-wave continuum with Fermi velocities v c and 



v s , respectively. Let M be the number of spin-down particles (er = —1) and N — M the 
number of spin-up particles for a total of iV on a ring of length L. We then twist jllj the 
Bethe-Ansatz equations obtained in Ref. |7j for N pseudo-momenta k = (k\, . . . , fcjv) and M 
rapidities A = (Ai, . . . , Xm) with a flux $ such that the total momentum P(k) = k = 0, 
i.e., 

M M N 

Lk = 2-Kl(k) - — $ + J2 0i{k - A) + J2 8o(k - k'), (3a) 

iV A k' 

M N 

= 2ttJ(A) + $ + #2(A - A') + - ( 3b ) 

A' k 

The two-body phase shifts for particle-particle, particle-spin wave and spin wave-spin wave 
scattering, i.e., 9 (k), 9i(k) and 9 2 (k) respectively, have been given previously M. The 
particle quantum numbers Ij and the spin-wave quantum numbers J a are integers or half- 
odd integers depending on the parities of N, M as well as on the particle statistics 0. 

In the thermodynamic limit, i.e., L — ► oo at fixed particle and spin wave densities, 
d = N/L and m = M/N, respectively, k and A will be distributed densely around the 
origin. The ground state is a filled Fermi sea characterized by the distribution function p(k) 
of particles and cr(A) of down-spins, so that we can replace the sums in Eq.(|3|) with integrals. 
This yields 

, 2-Klik) M a . fC n , f s f s , f B 

* ' /T\ i / (All/, . . \ — / . . \ w ■ ■ I / 



/ Bi(k - fi)a(n)dfji + f 9 (k - h)p(h)dh, (4a) 

27rJ(A) $ /-c r- 8 
= 7^ + t+/ 9 2 (\ - p)a(p)dp + 9 l (X-h)p(h)dh, (4b) 

and the values of S and C are fixed by the following equations: 

/ p{k)dk = d, (5) 

J — B 

j C a(X)d\ = md = (d-y)/2. (6) 

We rewrite these equations using the operator notation of ||: (k\Qi\h) = 9{(k — 
h), (k\p) = p(k), etc.. B and C are defined as projection operators to yield the finite 
limits of the integration. In this notation, Eq.((|) is 



k = + ®iCa + Q Bp, (7a) 

0=^^ + 1 + 8^ + 6^. (7b) 

Let us take the derivative of Eq.(|j) w.r.t. k. This then yields a set of integral equations 
for p{k) and cr(A) (Ref. |7[], Eq.(10)). In terms of the operator notation we have 

1/2tt = p + #iC<7 + K Bp, (8a) 



= a + X 2 C(T + K x Bp, (8b) 

where the new kernels are defined as (k\Ki\h) = w-O'^k — h). 

The neutral sector of the SC model is defined as M = N/2, i.e. C = II and y = 0. We 
now wish to calculate the effect of a small disturbence y ~ (C ^$ II) and $ ~ on the 
system. Eq.(8b") may be written as 



a + K 2 C(T = -K x Bp % (9) 
a = -(1 + J^Sp + (1 + J 2 )K 2 {1 - C)a, (10) 

with Jj the resolvent of -K^. Note that as pointed out in Ref. 0, the phase-shifts 9\ and #2 are 
essentially the phase shift of the H-I model after the identification of \s\ = Therefore 
the existence of the resolvents J\ and J 2 in the regime — 1 < s < has been proved in Ref. 

§■ 

Let o"o be a solution of Eq. fl8~15|) with C = IL, i.e., cr = —(1 + J 2 )KiBp. Then we may 
also write a + J 2 (l — C)a = (Jq. We substitute this back into Eq. fllaD and have 

k = — ^ + [0o - 6 X (1 + J 2 )Ki\ Bp - — $ - 0!(1 + J 2 )(l - C)(7. (11) 

Defining 8 = [6 - 8i(l + J 2 )lTi] and 5 = - 0^1 + J 2 )(l - C)a this yields 

fc = ^P + 05p + ^. (12) 



We see that we have succeed in rewriting Eq.(|3a|) as the equation at y = 0, <3> = plus a 
small perturbation term g. 

Let us briefly discuss how to handle the effect of this disturbance for additive quantities 
such as the energy. The net result of the presence of g will be a shift in the pseudo-momenta 
k = ho + 5k, with 8k of order 1/L. Then, expanding Eq.([T^), eliminating terms of order 1 
with the help of the unperturbed equation and defining 7 = 5kp, we have 7 + KB^f = gj2 r n. 
Introducing the resolvent (1 + J)B(1 + K) = (1 + K)B(1 + J) = II, we can solve for 7, i.e., 
7 = (1 + J)Bg/2ir. Let E = | YJ k 2 be the energy. Then the perturbed energy to first 
order will be E = YJ E(k + 8k) ~ YJ E(k ) + YJ E'(k )5k and we calculate the difference due 
to the presence of the perturbation g as 

AE/L = (E(k) - E(k ))/L ~ BE'Skp 
= [ B E'(kh(k)dk = ^ + BE' 

J-B 

= g+B(l + J)BE'/2ir = g + B(l + J)Bk/2n. (13) 

The excitation spectrum of Bethe Ansatz models has been expressed by Yang and Yang 
II | in terms of the solution to an integral equation. For the SC model the corresponding 



integral equation is 

e + KBe = k 2 /2 - (chemical potential). (14) 

The chemical potential is chosen such that e(±-B) = and thus we may take the derivative of 
the above equation, perform a partial integration and have e' + KBe' = k or e' = (1 + J)Bk. 
Inserting this into Eq. (|T3|) and performing another partial integration, we have the final 
result 

AE/L = -e + Bg'/2n = e+BK^l + J 2 )(l - C)a. (15) 

In Eq. (pb|) , J(A) was defined for A G [— C, C] only. However, the equations also work 
outside, i.e., = 2ttJ(oo) + $ + (M - l)0 2 (oo) + iV^(oo) = 2ttJ(oo) + $ + (M - 1)tt(1 + 
2s) - iV7r(l + s) ~ 2ttJ(oo) + $ + L[mdn(l + 2s) - dn(l + s)]. Thus we have 



a(X)dX = [J(oo) - J(Ci)]/L 

Ci 

~ [^(1 + s ) - ($/L)]/2tt, (16) 

since J(C\) = (M — l)/2 ~ Lmd. The same idea works for the lower limit. 

Let us now return to Eq.(|9]). We make the identifications a = ttX and p = ir\s\. Then 
the equation reads as 

2tx J-irC 2 cosh(a — p) — cos{2p) 2 J-b 2 cosh (a — irk) — cos(p) 

(17) 

If we compare this with Yang and Yang's || integral equation R + ~KR = £, we see that we 
can identify £ = —Ki(a)/2, K = /^(oO and write 

(7(A) = i / Bl i?(7rA - irk)p(k)dk. (18) 

2 J-B 2 

Furthermore, we have <7o = R^Bp/2 ~ RqBpq/2. Again from Yang and Yang, we can then 
compute the asymptotic behavior of ctq, i.e., 



<T — ► -e 2 M - / e 2 M p (k)dk 

U 2 J -Bo. 



IX _TTO_ 1 A 

-e 2m- / 

/i 2 j-s 2 

Ro{a)a/2, (19) 



where a is defined as the integral. Thus we write cr(A) ~ aR(irX)/2. Knowing the asymp- 
totics of R(a) from Yang and Yang M and using their notation, we can compute 



roo d roc 

/ a(X)dX = — / R(a)dat 
Jci 2n J-kCi 



= ^r(o), (20) 

with Ci = e -^ 2c i/ 2 ^ <^ i. The same calculation may be done for the lower limit and 
comparison with Eq.flTB]) yields 

Ci,2 = ^—W-p)T^/L)}. (21) 

Using Eq.fll^), we may then express the corrections to the energy as 



AE/L = -\t + BR Q {\ - C)a 



1 / e(k) \ I RohX - 7rk)a(\)d\ 

2 JB Ul-C 



dk 



-(TT 2 \-TT 2 k 



dk. 



(22) 



We now again use the identification a = tc\ and then shift the integration boundaries from 
Jc? U-lo) to / °° (/°oo). Furthermore, we treat the contributions from the upper and lower 
limits as independent. Then we are left with the Wiener-Hopf type problem of Yang and 
Yang || and we may conveniently use their results, i.e., 



AE/L 



" r e{k)e^ k '^dk- Und/fx) 2 + {^2/ v) 2 ] f (ztt/2//). 

7T L J 



4Jb 



(23) 



The k integration may be rewritten by a partial integration 

2 f i f jn.\^k/2u - 2 f i , 



- [ e(k)e^ k ^ = %\ e'(k)e^ k ^ = 
Jb n 2 Jb n 2 



(24) 



such that AE/L = £-T(iir/2fi) [( 2 + Q] -Using Eq.(|T]), we then have 



AE/L = t^Llltl h t 

' f 2 (0) ATT 3 



($/Lf + (tt - /i)V 



Again from Yang and Yang, we know T fl^ = s^iV-^) 



(25) 



and furthermore b/a = 2itv s . This 



gives the final result 



AE/L 



(tt - flf 



(26) 



Defining the stiffness D and the susceptibility x as AE/L = \ [x 1 (?//2) 2 + D{Q/L) 2 ] , we 
therefore have the desired relation 



(27) 



As we have shown in Ref. @, this relation may be also derived from conformal finite 
size formulas using the dressed charge matrix as calculated by Wiener-Hopf techniques [|7],|| . 
Furthermore, an analogous calculation yields the Luttinger relation for the H-I model and the 
method may be easily extended to other models. However, we made good use of the seminal 



papers of Yang and Yang in the present exposition and correspondingly useful results 

do not exist for most other models. We remark that the essential idea for our calculation is 
already contained in these papers, and also in Haldane's "Luttinger" papers |L|. 

In Ref. 0, we have also shown how to relate stiffness and susceptibility by threading 
the system with a flux of $ = 2tt(s + 1). The Bethe Ansatz equations @ for N particles 
and M down-spins simplify considerably at this point and are in fact the equations for iV 
particles and M — 1 down-spins at zero flux. This then gives D = x _1 /47r 2 (s + l) 2 . Thus 
we may express the susceptibility in terms of the spin wave velocity as \~ l = 27rv s (s + 1). 
In Fig.(|l|) we show a plot of x^ 1 as function of the interaction strength. The spin wave 
velocity has been calculated from the excitation spectrum of the SC model as given in Ref. 
||. As s — > the spin wave velocity approaches the velocity of a non-interacting single- 
component model, i.e. v s — > nd/2. Thus at half-filling, x _1 (0~) approaches the non-zero 
value 7r 2 /2 which is compatible with the result of Ref. |10|| . 

Finally, we wish to emphasize that we have derived the exact Luttinger liquid relation in 
the SC model starting from the asymptotic Bethe Ansatz equations. This further supports 
the use of the asymptotic Bethe Ansatz as an exact method in the thermodynamic limit. 

R.A.R. gratefully acknowledges partial support by the Germanistic Society of America. 
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FIGURES 

FIG. 1. The inverse susceptibility x" 1 for the SC model from the Luttinger relation for 
-1< s < 0. 
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